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■ CMS experiment is using Geant4 for Monte-Carlo simulation of the detector setup. Validation of physics processes 

f**^ ' describing hadronic showers is a major concern in view of getting a proper description of jets and missing energy for 

^Nj ■ signal and background events. This is done by carrying out an extensive studies with test beam using the prototypes 

^jq' or real detector modules of the CMS calorimeter. These data are matched with Geant4 predictions. Tuning of the 

, Geant4 models is carried out and steps to be used in reproducing detector signals are defined in view of measurements 

of energy response, energy resolution, transverse and longitudinal shower profiles for a variety of hadron beams over 
a broad energy spectrum between 2 to 300 GeV/c. 
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Series of Test-Beam measurements have been performed on the calorimetric system of CMS over the last few 
years in order to optimize its design and study its performance. Detailed Monte-Carlo simulations of the test-beam 
^ ' configuration have been made and the results compared with these test-beam measurements, leading to a number of 
• ^ improvements in the simulation code. Presented here is a comparison between the results from the 2006 Test-Beam 
and the Geant4-based Monte-Carlo simulations of this setup. 
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2. THE COMPACT MUON SOLENOID AND ITS CALORIMETRIC SYSTEM 

I The Compact Muon Solenoid detector (CMS)[l| is one of the general purpose detectors for the Large Hadron 
Collider (LHC) at CERN, to start operations in 2008. It's calorimetric system consists of Electromagnetic Calorimeter 
(ECAL) and Hadronic Calorimeter (HCAL) positioned inside the superconducting solenoid, and a pair of Very- 
QQ , Forward Calorimeters positioned outside of the return yoke of the magnet (figH^). For the purpose of this study only 
■ one super-module of ECAL crystals, without the preshower in front, and a pair of HCAL barrel wedges were used, as 
shown in figdjD. The Outer Calorimeter (HO) was present in both the test-beam environment and in the simulation 
• • i: geometry, but was intentionally left out of the comparisons, as its inclusion introduces additional uncertainties which 
. !_| ' were not well understood at the time of this study. The Forward Calorimeter (VF) was not included in the study. 

One of the two HB wedges used in the test-beam was wired in the normal "tower- wise" way (Fig 12^), while the second 
was wired in a "layer- wise" manner (FigOs) to allow for the studies of longitudinal shower profiles. Fig[3] shows the 
H2 beam-line arrangement at CERN's SPS accelerator. The VLE section enables the formation of very-low-energy 
beam tunes, allowing beam momenta as low as 2GeV/c to be studied. 



3. THE MONTE-CARLO SIMULATION 

The full simulation of the test-beam setup was performed using the CMSSW software frameworks, which uses 
internally the Geant4 toolkit Q , 3 • Fig|4] shows an overview of the detailed geometry used in the simulation. All 
beam-line elements were carefully described and their responses simulated. 

Several releases of the Geant4 code were used in the comparisons, and with each of them several physics-lists 
describing the interaction of hadrons with matter were tested. Historically the LHEP physics list was used as a 





Figure 2: Two different read-out scliemes for tfie HB wedges: Tower-wise (a) and Layer-wise (b) 



default, but it's limitations at high particle momentum warranted it's replacement with the QGSP physics list. 
Physics lists using the Bertini Cascade became available recently, bringing in new potential for better describing the 
interactions at low energies, and were used in this study as well. Most results in this paper are from release 9.1.p02 
of Geant4. 




400 GeV/c 
primary protons 



10 GsV/c < p < 400 GeVte 

secondary 
hadrons/electrons 



2-9 GeV/c lertiery 
hadrons/electrons 



Available beam tunes: 
pions 2-300 GeV 
muons 80/150 GeV 
electrons 9-100 GeV 



Wire Chambers (WC A,B,C): 
single hit to reject 
interaction in beam line 



P-ID: 

Cerenkov counter (CK2) - electron 

Cerenkov counter (CK3) - pion / kaon / proton 

Scintillators (V3, V6, VM) - muon tagging 



Figure 3: SPS's H2 beam-line. The VLE section is used for production/selection of very low beam energies (2 - 10 GeV) 
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Figure 4: Detailed MC simulation geometry 



9.0.p1 Response (MCideal calib.: eleSO) 
Birks': ECAL 0.1*BGO, HCAL 0.4*Def. 
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9.0.p1 Response (MCidealMIP calib.: eleSO) 
Birks': ECAL O.I'BGO, HCAL 0.4*Def. 
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Resolution (MCidealMIP) 
Birks': ECAL O.fBGO, HCAL 0.4*D6f. 
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Figure 5: Response (left) and resolution (right) of the simulated calorimertic system compared to the measured in TB2006 
for the case of combined ECAL+HCAL (top) and HCAL alone (bottom). 



4. RESULTS 

4.1. Linearity of response and Resolution 

FigO shows the comparison of linearity of response and resolution for the combined system ECAL+HCAL, as 
well as for the HCAL alone. Overall the agreement is good, except for resolution of the simulated detector being 
significantly better at high energies. The source of that discrepancy can be identified in the significantly narrower 
distribution of HCAL signal at high energies in FiglHl However, the general agreement of the simulated spectra and 
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ECAL+HCAL (left); ECAL (center); HCAL (right); 50GeV (top); and ECAL+HCAL system for 2GeV (a); 20GeV (b); 
200GeV (bottom) 50GeV (c); and 200GeV (d); pion beams 

Figure 6: Reconstructed energy spectra of simulated events compared to TB2006 data 
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Figure 7: Longitudinal shower profiles for SOOGeV (a); lOOGeV (b); and 30GeV (c) pion beams compared to data. 



the measured ones is clear. 



4.2. Shower Profiles 

FigEI shows the comparison of the longitudinal shower profiles measured in HB2 wedge and simulated with three 
different physics lists. The QGSP-BERT physics list, using the Bertini-Cascade for low-energy interactions, shows 
best agreement with the measured profiles. Fig[5] compares the transverse shower profiles measured in ECAL, and 
the transverse shower containment. Agreement is good for both physics lists shown. 



4.3. Fraction of MIP-lil(e events 

An important measurement for the calibration of HCAL is the fraction of events which leave MlP-like signal in 
ECAL. As Figini shows, the Bertini physics list has a major discrepancy at low energies, which is attributed to the 
inadequate description of the quasi-elastic processes. This is one of the major issues which the Geant4 team is 
addressing in the next releases of the software. 
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Figure 8: Comparison of ECAL Transverse shower profiles in 77-direction (a) and (^-direction (b) for lOOGeV pions tliat 
interacted hadronically in the crystals; and shower containment of these showers (c) as a function of crystal matrix size (NxN), 
normalized to beam energy. 
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Figure 9: Fraction of pion events leaving MlP-like signal in ECAL, as a function of beam energy, compared to data. Super- 
imposed is the same quantity for muon beam, and muon beam scaled down to match pion data at SOOGeV 

5. CONCLUSIONS 

Geant4 simulated response is in overall agreement with the measured data. We still observe small discrepancies in 
the following quantities: 

• Simulated resolution is much better than the measured one, which is due to the smaller variation in signal 
distribution in HCAL at high energies. 

• Some discontinuities in the response curve at intermediate energies. 

• Fraction of MlP-like events in ECAL has peculiarities due to quasi-elastic processes in Bcrtini models. 
We continue our active collaboration with the Geant4 team to fully understand and correct these discrepancies. 
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